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The self-consistent determination of the damping rate of a fast moving fermion in a hot QED 

I plasma is reexamined. We argue how a detailed investigation of the analytic properties of the 

£f) . retarded fermion Green's function motivated by the cutting rules at finite temperature may resolve 

' ambiguities related to the proper definition of the mass-shell condition. 

m 
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q^. I. INTRODUCTION 

<D ' 

Recent studies of damping rates of fast moving particles in a hot QED (or QCD) plasma have raised some interesting 
problems [1-11]. The difficulties start from the infrared sensitive behaviour of these rates: a logarithmic divergence 
remains even after including the perturbative Braaten-Pisarski resummation 12[l3|] , which only screens the infrared 



sensitivity down to scales of O(eT). In the case of QED, which has no magnetic mass to serve as infrared cut-off at 
the scale 0(e 2 T), a self-consistent determination of the rate 7 of a fast moving fermion in the plasma is the most 
elegant solution, as first suggested by Lebedev and Smilga Q . Here 7 also plays the role of the infrared cut-off. In the 
analysis considered until now it is always assumed (even explicitly) that the retarded propagator of the fast fermion 
(more precisely its analytically continued form) has a true complex pole in the lower energy half plane at the position 
E — £7, with E the energy of the fermion. However, the presence of this pole does not allow an infrared stable solution 
for 7 when the on-shell condition at this complex pole is required ||: this is indeed the favorable and consequent 
condition in the case of a true pole on the "physical" Riemann sheet! 

Based on the analytic structure of the retarded Green's functions as deduced from the general properties of spectral 
functions we argue in the following that the complex pole in question is actually not on the "physical" sheet, and 
we construct an explicit and simple example which should represent the physical (realistic) situation and which may 
help to clarify the problem of self-consistency for 7 [Q . By this attempt we take the point of view that the damping 
rate of a fast moving fermion in a QED heat bath is a physical quantity, contrary to the arguments given in Ref. [|| . 
Therefore an infrared finite result is to be aimed for it. 
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II. RETARDED GREEN'S FUNCTION AT FINITE TEMPERATURE 



First it is convenient to illustrate some properties of the retarded Green's function G R (po,p) [16-20]. It is defined 
as the boundary value of a complex function G(z,p) as the complex argument z approaches the real axis from above: 



G R (Po,p) = G(z = po + ie,p), 
where the function G{z,p) is determined by a spectral function p{po,p) as 



G(z,p) 



dp'p 
Po-' 



P(Po>P)- 



If p{po,p) is a Lorentzian of the form 



p(Po,p) 



7/tt 



7 >0, 



( Po -E(p)r+^-' 

for all values of po, — oo < po < +00, then the retarded Green's function becomes 

iRu ^_ f°° d Po 7A 



(1) 



(2) 



(3) 



G h (po,p) 



p' -(po + ie) (p' - E{p)f + 7 2 
1 



Po + ie- (E(p) - ij) ' 



(4) 



where we first keep po on the real axis. However, this explicit function may be analytically continued to complex 
values of po : for Im po > G R is regular as it should be from the definition (1), but for Im po < the analytically 
continued G R has a true pole at po = E — i-f. The origin of this pole is easily understood when considering the 
situation in the pp-plane 1 21 Hi- As illustrated in Fig. la for the continuous integration contour — 00 < p' < +00, 
the singularities of the integrand are at p' = po + K and due to the Lorentzian at p' Q = E ± i-f, indicated by crosses. 
The contour for the retarded function is below the pole at p' = p + ie, indicated by e > 0. In order to continue G R 
into the lower half plane, Im po < 0, the contour has to be deformed (Fig. lb), and when p' = Po + ie approaches 
p' Q = E — i-f it is pinched, resulting in the pole of G R at po = E — i-f. For the continuous contour (Fig. 1) this pinch 
is always present when the continuation to values for Im pq < is performed. 

However, the situation becomes different when, for instance, it is assumed that the spectral function p is non- 
vanishing at some threshold po > p t h, but otherwise it is zero. If, for example, p is given for po > p t h by Eq. (^|) with 
E(p) > Pth, the retarded function is 



G R ( Po ,p) 



dp' 



7/tt 



Pth Po 

JLf. 



(Po 



(P'o 



2iri I po - (E(p) - ij) 



In 



E{p)f+ 1 2 



Pth -po-ie 



po - (E(p) + i-f) 



In 



Pth - E{p) + ij 

Pth -po-it 
Pth - E(p) - i-f 



(5) 



In the Po"Pl ane the situation looks as in Fig. 2. One immediately sees that when p'o = po + ie approaches p' = E — i-f, 
for example, the pg-integration contour is not necessarily pinched, and so there is no pole at po = E — i-f in the 
analytically continued G R on the first ("physical") sheet. Instead, there is an endpoint singularity at po = pth, which 
is a logarithmic branch point, and the discontinuity across this cut on the real axis is given (by construction) by the 
Lorentzian p{po,p) of Eq. (||). With respect to this cut the pole has been moved onto the other sheets obtained after 
continuation of the logarithm, ln(p t h — Po — ie): this amounts to deforming the contour in such a way that a pinch is 
present (cf. Fig. lb). 

The behaviour with respect to po — E + i-f is symmetric with respect to the one at po — E — i-f; on the first sheet 
G R (po,p) of Eq. (^|) is regular for Im po > and for Im po < 0, with the only singularity on this sheet that of the 
logarithmic branch point at po = pth- 

The case just described is in close analogy to the discussion of resonances (at zero temperature): with respect to 
the variable of the energy squared resonances are on the unphysical sheet, and the branch point is determined by the 
threshold properties of scattering amplitudes. 
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III. FERMION DAMPING RATE 



As the simplest case we consider the damping rate of a heavy fermion ( "muon" ) of mass M in a hot QED plasma 
of (massless) electrons and photons. The energy, the momentum and the velocity of the "muon" are denoted by 
E,p (p = \p\), and v, respectively, but the main interest is in the limit v — ► 1 |B,B[LC3 . 

Since we are only concerned with the leading order behaviour e — > at high temperature T, for a fast muon 
(E 3> M > T) a couple of approximations which simplify the calculations and the discussion may be applied [1-10]. 
With -{G R )- X =p - S fl and Eq. (g) we relate in the usual way 7 to the imaginary part of the muon self-energy by 
l(Po,p) = -ImS fl (po,p), with p real. The explicit expressions are derived, for example, in detail in jn|, which we 
closely follow concerning conventions and notation, and we start off with Eq. (17) of this reference. To leading order 
mlwe have 

l(Pt),p) - +e 2 T [ 7I-T3 / — pt(qo,q)lmG R {p - qo,p- q). (6) 

J soft ( 27T ) J-q 90 

The following remarks summarize the results for 7 obtained so far in the literature: 

- Eq. (|^) is derived in the one-loop approximation, in which the hard energetic "muon" emits/absorbs one soft 
photon of four momentum g M . As well, corrections of order 1/E are neglected. 

Only the dominant transverse photon contribution has to be taken into account, which following the Braaten- 
Pisarski hard thermal loop resummation method jl2| is determined by its spectral density pt (qo , q) ; its explicit 
form is given in p3| , but for the following we only require the (approximate) integral 

[ +q ^p t (q ,q)^l/q 2 (7) 

J-q qo 



when q < eT §,|JT0|. 

All the complications due to the spin of the fermion are suppressed in Eq. 
fermion propagator is described by a (retarded) scalar function G R . 



), in the sense that the heavy 



Inserting for a bare fermion —(G R ) 1 = po + ie — E(p), E(p) = \fp 2 + A/ 2 , into Eq. (||), we have 

lmG R (p - qo,p- q) 



90—0 



nS( Po - v/(p-g) 2 +M 2 ), 



(8) 



and one immediately finds the infrared divergent result 

T I 

4ir 



j(po — E,p) ~ ^—T I (9) 



This is mainly due to soft photon exchange, and it is the origin of the problems with the hard fermion damping 
rate. It also shows that to leading order only the infrared sensitive region |go| < q — * has to be considered. It 
is worth noting that the coefficient in Eq. (g) is gauge parameter independent [||j|] . 

- In hot QED there is no (non-perturbative) magnetic mass to provide a cutoff to the logarithmic infrared diver- 
gence in Eq. (||) [Q. In QCD the magnetic mass is expected to be on the scale m mag ~ g 2 T, with g the strong 
coupling constant: consequently using m ma „ and evaluating 7 on the real axis (po — E) a finite value may be - 
and has been - argued for quarks in QCD ||,[4],f7j . 

- From Eq. ([)]) 7 is "anomalous" , in that its magnitude is on the scale e 2 T (neglecting lne factors for the moment). 
Therefore it has been first conjectured by Lebedev and Smilga || to use 7 itself as a possible infrared cutoff. 
This opens the possibility of a self-consistent calculation of 7 by replacing G R in Eq. (||) by the simple Lorentzian 
of Eq. (||) j|. In Ref. Q it is shown that only the fermion propagator should be modified, but not the photon 
one; also no vertex corrections are required. 
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- The self-consistency requirement for the damping rate requires, however, clarification of the "on-shell" condition 
used to evaluate Eq. @. On the one hand, as used in Ref. [g], one can keep in Eq. ([|) p on the real axis, in 
which case po = E(p) is used. On the other hand, one can generalize Eq. (||) to complex po by demanding that 
it holds at the complex pole po — E{p) — under the narrow width assumption 7 <C E{p). It is crucial to point 
out that under the assumption that G R is given by Eq. ([|) the complex pole is on the first (and only) sheet 
in the energy plane; therefore the second case of self-consistency could be argued as the physical one, and is in 
any case a point at which gauge invariance can formally be proven to hold [ B4J . However, as remarked in Ref. 
||, in this case the infrared divergence is not screened by a non- vanishing 7! Therefore, this attempt fails for 
QED when Eq. (@) is used as a model for the "dissipative" retarded fermion propagator, and as such the narrow 
width condition in the form Im£^(po — E) ~ ImE fl (po = E — 17) does not hold. 

Rather than concluding that the fermion damping in QED is an unphysical quantity, and therefore not observable 
as argued in Ref. Q , we take the point of view that the ansatz of Eq. (Q) for the retarded energetic fermion propagator 
G R does not reflect the proper physical conditions. Instead it is realistic to use as a simple model of the spectral 
function the form given by Eq. (||), which allows for branch cuts on the real axis in the energy plane. 



IV. SPECTRAL FUNCTION AND CUTTING RULES 



When discussing the location of branch points and cuts in self-energy functions at nonvanishing temperatures it is 
useful to recall that Weldon ]25[ | has shown, starting from the one-loop approximation and extrapolating to the many 
particle case, that the branch points are determined by the T = masses of the particles in the heat bath. One can 
see this by considering the spectral function directly, where one notes that it is evaluated - independent of T = or 
T ^ - by using the energy eigenstates of the (full) Hamiltonian, H\n >= E n \n >. For example, for fermions the 
resulting p at finite T has the following structure ]Tfy|T8| ]: 

P(P0,P )==X^ 27r ) 4<5 ( P ° ~~ ( En ~ E m))3(p- (Pn -Prn)) 
m,n 

. e - E ™/ T (l - e- p °/ T )\(m\Tjj\n)\ 2 , (10) 

where ip denotes the Dirac field operator. Obviously the energies E n and momenta p n appearing in the (5-functions 
arc temperature independent, and therefore so are the positions of the branch points according to the cutting rules 

HI!- 

Although the positions of the cuts do not depend on temperature, the discontinuities across the cuts become 
temperature dependent. Let us consider the one-loop g'ipiprf' self-energy example of Weldon [ p5[ . The cut structure 
for the fermion self-energy with mass M is reproduced in Fig. 3 in terms of the variable s — p\ — p 2 - as at zero 
temperature a small (zero temperature) photon mass A is introduced, having in mind the limit A — > whenever allowed. 
The cut starting at s > (M + A) 2 is familiar from zero temperature. The cut between — (M 2 — A 2 ) < s < (M — A) 2 
is due to the absorption/emission of photons from the heat bath, and vanishes for T — > 0; in the limit M > T its 
discontinuity is exponentially suppressed, and therefore this cut is neglected for the following discussion. In Fig. 4 
we plot the discontinuity in this example for both hard and soft regions of external momenta - we consider the two 
contributions ^4(po;P)7oPo and MD(po,p) to ImI] R (po,p). We have kept the photon mass A small but finite in these 
figures in order to differentiate between the two regions s > [M + A) 2 and s < (M — A) 2 - this is indicated by the 
break in the curves. In the soft regime of Fig. 4(a) we find that the Landau damping contribution for ,s < (M — A) 2 
dominates, as expected, but in the hard regime of Fig. 4(b) the cut coming from s > (M + A) 2 starts to dominate. 

From this we deduce that the thermal spectral function for the heavy, energetic fermion propagator — (G^)" 1 = 
Po — ^ R (po,p) has to have a contribution from at least the branch cut for s > (M + X) 2 , i.e. for \p \ > \Jp 2 + (M + A) 2 . 
Note that, in this region, such a contribution comes entirely from the self-energy H R , since there is no pole contribution. 
Consequently we take as a realistic ansatz for the determination of G R in Eq. (0) the following: 



(i) for positive energy there is a single cut starting at pth = \Jp 2 + (M + A) 2 ; this is smaller than the energy of 



the thermally excited heavy fermion, which receives contributions of O(eT) fl23| , p7 28 1, such that E(jp) — pth 
0(e 2 T 2 /E) for A — > and E > T. Except for this, in the following the thermal contribution of O(eT) to the 
mass is neglected because of the limit M > T; 

(ii) because of the narrow width condition 7 <C E(p), near po ~ E(p) the cut's discontinuity is dominated by the 
nearby "pole" . We assume the pragmatic parametrization of this discontinuity is given by the Lorentzian of 
Eq. (|3|), assuming 7 to be momentum independent (as a calculational simplification); 
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(iii) in order to respect the symmetry properties of the fermion spectral density p9[ , p{—po,p) 
following ansatz is suggested: 



p(Po,P) 



7 



2ttE( P ) 



1 



1 



(p -E(p)) 2 + 7 2 (p + £(p)) 2 +7 2 



0(bo| -pth)- 



-p(Po,P), the 



(11) 



The temperature dependence only shows up in 7 = j(T); 

(iv) under the strong assumption - which we accept for simplicity in the following - that Eq. ([ll]) dominates p for 
all values of po, and not only in the neighbourhood of E(p), one may require the sum rule to be satisfied by the 
ansatz (pi]): 



r+00 

1 = / p(po,p)p dp ; 



(12) 



this fixes the normalization given in Eq. ([Ll]) when terms of O^/E) are neglected. This assumption leads to 
an overestimate of the damping rate when determined self consistently 

Considering these points, we take as the retarded Green's function for this simple toy model that following from the 
spectral function of Eq. (|ll|) M: 



G R (po,p) 



-(7 -» -7) 



1 



AnE(p) 



Pth + E(p) - ij p th - p - it 
Po + ie- E(p) + ij \p th - E(p) + ij p th + p + ie 



■In 



E{p) -» -E{p) 



(13) 



This function - more precisely its analytic continuation in po - does not have poles at po = ±E ± 17 on the first 
"physical" sheet with respect to the branch points at po = ±pth- The discontinuity across the cut starting at 
Po = ±p t h is given - by construction - by the Lorentzian of Eq. ([ll]) for real values of po, and near po — E{p) we find, 
for 7 <C E(p) and E > T, the narrow width condition 



lmG R (po = E(p)) 
lmG R (p = E(p)-ij) 



1. 



(14) 



V. SELF-CONSISTENCY FORMULATION 



Although we have assumed the spectral function of the propagator has the form of Eq. ( |1 1| ) for all values of po, in the 
true situation we might expect that this would be only in a neighbourhood of po — E{p). Thus, for a self-consistent 
determination of 7 we insert G R of Eq. (|l^) into Eq. (||) and evaluate the result at the point po — E(p); to leading 
order we then have 



e 2 f eT f +1 
7(p - E(p),p) ~ J dq J dcosf 



7 



eT r + l 

rl nr\G f, 

7 2 + q 2 cos 9 

e 2 f eT dq , , . e 2 f eT dq 

—T — arctan(g/7) ~ — T / 
2tt 2 Jo q 4tt J 1 q 



e 2 eT e 2 1 
~^Tln— ~— Tin- (15) 

47T 7 47T e 

where the infrared "screening" by 7 is explicitly exhibited. This then reproduces the original self-consistent derivation 
of the fast damping rate in Ref. [g| , but in this case without a singularity if one had used instead the point po = E — ij 
on the physical sheet to evaluate Eq. (|^) . This is consequently consistent with the narrow width assumption of the 
form lmT, R (p = E(p)) ~ lmY, R (p Q = E(p) - fy). 

With minor modifications the preceding mechanism should be applicable to the case of QCD fast damping rates 
and colour relaxation times [pl[ , without having to introduce a magnetic mass as an infrared cut-off. 
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VI. LARGE TIME BEHAVIOUR AND DISCUSSION 



In order to interpret 7 of Eq. (|l5|) in the context of this toy model, we study the time dependence of the Green 
function G R (t). This is given by the Fourier transform of the retarded Green's function of Eq. (|l3|), 

G R (t) = / d Po G R ( Pa )e-^ t , (16) 



where we now suppress the reference to the spatial momentum p. Using the spectral representation of G R (po) we 
obtain for times t real and positive 



/ + 00 p + 00 —ip t 

r+oo 

= 2m / dp Q p(p Q )e-^ t . (17) 



Inserting now the spectral function under consideration, Eq. (pd|), into Eq. (|T^), we find 
G R {t) = 2m-L- \ dp - — ^—3 + c.c. 



2wEj pth ^( Pa -E{p)f+^ 
— { Trsm{Et)e-^ - / dx- LV 



E I J(E- Pth )h x2 + 1 J 

— {tt sin(M)e" 7 * + cos(^) [e" 7 *£;i(7i) - e 7t K(- 7 t)] } , (18) 

where Ei{x) is the exponential integral JB2| and the approximation E 3> T, for which (E—p t h) 7 and [E+pth] 3> 7, 
has been used. One can consider Eq. ( |18| ) in two limits - if we assume jt S> 1 then we find 

G , (t) .|^) +0 (1) 2 , 

while if we assume 7 t <C 1 then we obtain 

G R (t) - ^sin(^i) + ^ [2(ln( 7 <) + 1 - j E ) cos(Et) - irsin(Et)} + 0(~/t) 2 , (20) 
E E 

where je is Euler's constant. Thus, as noted in Rcf. JllJ and stressed in Ref. g, the time dependence of G R (t) 
may not necessarily be of an exponential form, even for very large times, and so care must be taken in these cases in 
characterizing 7 as an exponential "damping" rate. Even so, it is still a parameter within the context of the ansatz 
for the Green function which remains to be determined, and for this we can use the self-consistent condition derived 
from Eq. (|Tl|): 7 ~ — ImS R (po — E); the relation ImY, R (po — E) ~ ImE fl (po = E — 17) assures us to this order that 
using the "complex" on-shell condition po — E — i 7 will lead to the same self-consistent determination of 7. 
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FIG. 1. Integration contour for the retarded Green's function. Poles of the Lorentzian form of Eq. (|3|) are indicated by (g>. 
In (a) po is real, while (b) shows the pole singularity for Im po < which arises from the pinching of the contour between the 
two singularities p' = po and p'o = E — i*y of the integrand. 

FIG. 2. Integration contour for Eq. (jEJ), p t h < Po < oo, indicating that there is no pinch singularity on the first sheet at 
Po = Po — E =p ij with respect to the endpoint singularity at po = pth- 

FIG. 3. Location of the branch cuts (in the one- loop approximation) for the thermal fermion propagator EfJ . A denotes the 
T = photon mass (A -» 0), while M is the fermion mass. 



FIG. 4. The discontinuity of the one-loop fermion self-energy example of Weldon j25|. In (a) we consider the soft region, 
with M/2T = 0.2 and p/2T = 0.1, while in (b) we consider the hard region, with M/2T = 15 and p/2T = 10. The upper line 
in both figures denotes A(po)/A(po = oo), while the lower line denotes D(po)/ D(po = oo). 
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